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Many single molecule experiments for molecular motors comprise not only the motor but also 
large probe particles coupled to it. The theoretical analysis of these assays, however, often takes 
into account only the degrees of freedom representing the motor. We present a coarse-graining 
method that maps a model comprising two coupled degrees of freedom which represent motor and 
probe particle to such an effective one-particle model by eliminating the dynamics of the probe 
particle in a thermodynamically and dynamically consistent way. The coarse-grained rates obey 
a local detailed balance condition and reproduce the net currents. Moreover, the average entropy 
production as well as the thermodynamic efficiency is invariant under this coarse-graining procedure. 

Our analysis reveals that only by assuming unrealistically fast probe particles, the coarse-grained 
transition rates coincide with the transition rates of the traditionally used one-particle motor models. 
Additionally, we find that for nrulticyclic motors the stall force can depend on the probe size. We 
apply this coarse-graining method to specific case studies of the Fi-ATPase and the kinesin motor. 


I. INTRODUCTION 

In many single molecule experiments beads that are at¬ 
tached to molecular motors are used to infer properties 
of the motor protein from the analysis of the trajectory 
of these probe particles. In particular, external forces 
can be exerted on the motor via such a probe particle 
[l], 2j. In the theoretical analysis of such assays, the mo¬ 
tor is usually modelled as a particle hopping on a discrete 
state space with transitions governed by a master equa¬ 
tion [SHa]. Alternatively, the so called ratchet models 
combine continuous diffusive spatial motion with stochas¬ 
tic switching between different potentials corresponding 
to different chemical states [1, [bij • These approaches of¬ 
ten comprise only one particle explicitly, representing the 
motor. The contribution of external forces which in the 
experiments act on the motor only via the probe are then 
included in the transition rates @,[6, 11 17] (or Langevin 
equation for the spatial coordinate [lg, L!jl ) of the motor 
particle directly. However, theoretical models that are 
used to reproduce the experimental observations should 
comprise at least two (coupled) degrees of freedom, one 
for the motor and one for the probe particle. Such models 
consisting of one degree of freedom hopping on a discrete 
state space representing the motor coupled to a contin¬ 
uously moving degree of freedom representing the probe 
are discussed in |2I> |2?|. While multi-particle models are 
more precise and better represent the actual experimen¬ 
tal setup, one-particle models are widely-used toy models 
often applied to illustrate basic ideas. 

Simplifying the description of systems consisting of 
many degrees of freedom with a concomitant large state 
space while still maintaining important properties is com¬ 
monly known as coarse-graining. In the context of 
stochastic thermodynamics (2j|, various coarse-graining 
methods have been applied, e.g., lumping together states 
of a discrete state space among which transitions are fast 
averaging over states for discrete [H] or contin¬ 
uous processes |32|, |34j , eliminating single states from a 


network description 1351 - 1371 or eliminating slow (invisi¬ 
ble) degrees of freedom [38l - t40| . It was found that, in 
general, coarse-graining has implications on the entropy 
production and, in particular [41| . dissipation. In the 
context of biological systems and especially molecular 
motors, coarse-graining procedures mostly focus on elim¬ 
inating selected states of the motor [37], j42| or on reduc¬ 
ing continuous (ratchet) models to discrete-state models 

mnn. 


In the present paper, we introduce a coarse-graining 
procedure that allows to reduce molecular motor-bead 
models to effective one-particle models with discrete mo¬ 
tor states with the external force acting directly on the 
effective motor particle. We eliminate the explicit dy¬ 
namics of the probe particle completely still maintaining 
the correct local detailed balance condition for the ef¬ 
fective motor transition rates and preserving the average 
currents of the system. As a main result, we find that the 
coarse-grained rates show a more complex force depen¬ 
dence than the usually assumed exponential behaviour 
and a more complex concentration dependence than mass 
action law kinetics. 


The paper is organized as follows. In section El we 
introduce our coarse-graining method on the basis of a 
simple motor-bead model with only one motor state and 
apply it to a model for the Fi-ATPase (26| . In section UTTl 
we generalize the procedure to motor models with several 
internal states and apply it to both a refined model for 
the Fi-ATPase and to a kinesin model. A possible experi¬ 
mental implementation of our method is presented in sec¬ 
tion ED We show that entropy production and efficiency 
remain invariant under this coarse-graining procedure in 
section 0 discuss implications on the stall conditions in 
section El and conclude in section Em 
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II. GENERAL ONE-STATE MOTOR MODEL 
A. Explicit motor-bead dynamics 

The general model for motor proteins with only one 
chemical state consists of one degree of freedom repre¬ 
senting the motor which jumps between discrete states 
n{t) separated by a distance d. The motor is coupled 
with the second degree of freedom representing the probe 
particle via some kind of elastic linker, see Fig. [T] [26j . 
The motion of the probe particle with continuous co¬ 
ordinate x(t) is described by an overdamped Langevin 
equation with friction coefficient 7 and constant external 
force / ex , 

x{t) = (~d x V(n -x) - / ex )/7 + C(*), (1) 

including the potential energy of the linker V(n — x) 
and thermal noise C,(t) with correlations (C(^ 2 )C(^i)} = 
2<5(i 2 — t\)h- Throughout the paper, we set k^,T = 1. 
This choice implies that the product of force / ex and 
distance d appearing in the figures below is measured 
in units of k&T. The (instantaneous) distance between 
motor and probe is denoted by y. The system is charac¬ 
terized by the pair of variables ( n,x ) and is “bipartite” in 
these variables since transitions do not happen in both 
variables at the same time. The transition rates of the 
motor fulfill a local detailed balance (LDB) condition 

W , ^ = exp[A/i - V(y + d) + V fa)]. ( 2 ) 

w (y + d) 

The free energy change of the solvent A/j = /it — pn — AT 
with )ii = Hi q + ln(cj/c® q ) and nucleotide concentrations 
Ci is associated with ATP turnover. The probability den¬ 
sity p(y) for the distance y obeys a Fokker-Planck-type 
equation 

dtp(y ) = d v {(d y V{y) - / ex ) p(y) + d y p{y)) /y 

+ w + (y - d)p(y - d) + w~ (y + d) p(y + d) 

- (w + {y)+ w~(y)) p(y). (3) 


w (n,x) w + (n,x ) 



x(t) n(t) 

y(t ) 


FIG. 1. (Color online) Schematic representation of a motor- 
bead model comprising a one-state motor (small blue sphere) 
attached via an elastic linker to the probe particle (large red 
sphere). An external force / ex is applied to the bead. The 
transition rates of the motor are denoted by w + ( n , x ) and 
w~(n,x). The load sharing factors 9 + and 9~ indicate the 
position of an underlying unresolved potential barrier relative 
to the minimum of the free energy landscape of the motor. 

particle hopping between states separated by d. We 
thus have to eliminate the x-coordinate from the ( n,x)~ 
description resulting in a system characterized only by 

n. 

For the coarse-grained model, we impose the following 
conditions. The coarse-grained transition rates which 
advance the effective particle by d should obey a LDB 
condition 

57+ 

— = exp[Ap, - fexd] (5) 

as the force is now assumed to act directly on the effective 
motor particle. Furthermore, we require that the coarse¬ 
grained particle moves with the same average velocity 
in the steady state as the motor and the probe in the 
original model, i.e., 

v = d(0+- FT). (6) 

Solving the linear system of equations ( 0 © yields the 
coarse-grained rates 


For constant nucleotide concentrations, the system 
reaches a non-equilibrium stationary state (NESS) with 
constant average velocity 


v = d p s (y)(w + (y) - w (y))dy 
J — OO 

/ OO 

P s (y) ( d v V(y ) /ex) /7 dy 

-00 

and stationary distribution p s ( y ). 


( 4 ) 


B. Coarse-graining procedure 

In the coarse-grained description of the model we want 
to map the motor-bead system to one effective motor 


n+ = 


rr = 


vexp[Ap - f ex d]/d 
exp[A/z - / ex d] - 1 
v/d 

exp[A/z - / ex d] - 1' 


( 7 ) 

( 8 ) 


The coarse-grained rates can be interpreted as effective 
transition rates that correspond to a transition process 
after which both particles, motor and probe, have ad¬ 
vanced a distance ±d. In principle, there are (for any y) 
many possible displacement processes to advance both 
particles by d, including ones with l forward and l — 1 
backward motor jumps. The coarse-grained rate corre¬ 
sponds to the rate with which one such effective displace¬ 
ment will happen. 

In general, the coarse-grained rates depend (via v ) on 
all model parameters, including the friction coefficient of 
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the probe particle and the specific potential of the linker. 
If one had chosen coarse-grained rates by just averaging 
over the positions of the probe particle, i.e., by 

/ OO 

p s (y)w ± {y)dy, (9) 

-OO 

one would have obtained rates that yield the correct av¬ 
erage velocity but do not fulfill the LDB condition, as 
discussed in section HlEl below. 

For a more explicit analysis, we must specify the for¬ 
ward and backward rates of the motor. We choose j26j 

w + (y) = wo exp[/i + - V(y + d9 + ) + V(y)\ (10) 

w~(y) = wq exp[/i“ -V(y- dd~) + V(y)] (11) 

where 9 + and 9~ are the load-sharing factors with 
9 + + 9~ = 1 and y + = y~ = yp> + yp- We assume 
an exponential dependence of the transition rates on the 
potential difference of the linker according to Kramers’ 
theory. This exponential dependence on the potential dif¬ 
ference is similar to one-particle models where the rates 
of the motor typically depend exponentially on the exter¬ 
nal force with a corresponding load-sharing factor si. 

C. Time-scale separation 

In this section, we will investigate under which condi¬ 
tions the coarse-grained rates (0® can be expressed us¬ 
ing a single exponential dependence on the external force 
as typically assumed for mechanical transitions within 
one-particle models Si- 

Inserting eqs. mm in eq. © in the NESS shows 
that the contribution due to motor jumps is weighted 
with a (dimensionless) prefactor 

e = wq exp[/r^ q ]d 2 7 - ( 12 ) 

Here, wo exp[/x!j?] determines the timescale of the transi¬ 
tions of the motor while yd 2 determines the timescale of 
the dynamics of the probe particle. The latter is mainly 
governed by the size of the bead and the step size of the 
motor whereas Wq exp [//()?] is determined by the attempt 
frequency and also by the absolute nucleotide concentra¬ 
tions. 

If the dynamics of the bead is much faster than the 
transitions of the motor, time-scale separation holds with 
e -> 0 |3ll, S3- In this limit of fast bead relaxation, 
denoted throughout by a caret, the stationary solution 
of eq. © in the NESS becomes 

p s (y) = exp [~V(y) + f ex y]/J\f (13) 

with J\f = f exp[— V(y) + / ex y] d y. The average veloc¬ 
ity is then given by 

/ OO 

p s {y){w + (y) -w~(y))dy 

-OO 

=dwo (> T - / “ d0+ - e» D+f * p+f °* d0 ~') . (14) 


This expression inserted into eqs. ©0 yields 

Ci+= Wo e llT - f ^ de+ , (15) 

fr = w 0 e liD+ltp+f ” M ~ (16) 

independent of any specific linker potential V(y). Since 
this force dependence is purely exponential with the cor¬ 
rect load sharing factor, these expressions represent ex¬ 
actly the rates typically used in one-particle models. We 
notice that within this approximation = ( w + (y )) and 
= (w~ (y)) holds true, which is in agreement with 
other coarse-grai ning pr ocedures in the time-scale sepa¬ 
ration limit, e.g., [3lM33l | . 

Note that only transition rates of the motor whose de¬ 
pendence on the linker potential is chosen accordingly in 
the Kramers form (eqs. (fTOlfTTl) ') lead generically to con¬ 
sistent coarse-grained and averaged rates when using the 
fast-bead limit of p s (y). 

D. Example: Fi-ATPase 

In general, a strong time-scale separation between mo¬ 
tor and probe is not necessarily realistic. In this case, eq. 
© must be solved numerically. We will use the model 
introduced in [2(1 . see Fig. © with a harmonic poten¬ 
tial V{y) = ny 2 /2 as a simple example to illustrate our 
coarse-graining procedure. 

In Fig. © the results for fl + and f l~ are shown for vari¬ 
ous values of the friction coefficient 7 . With decreasing 7 , 
the rates approach their corresponding fast-bead limits, 
f2 + and fl - . These values are upper bounds because de¬ 
creasing 7 implies smaller probe particles which exert less 
drag on the motor. For finite 7 , the coarse-grained rates 
do not show a single exponential dependence on / ex over 
the whole range of external forces. Such a dependence, 
however, is usually assumed to hold within one-particle 
models. Moreover, the coarse-grained rates depend on 7 , 
which is a parameter not incorporated explicitly in many 
one-particle models. 

The experimentally accessible values of 7 cover a wide 
range of the values chosen in Fig. ® A dimer of 
polystyrene beads (~ 280 nm) as used in (Tol (52| cor¬ 
responds to 7 = 0.5 s /d 2 (red (dark gray) line with 
squares) while a 40 nm-gold particle |52l - l54j corresponds 
to 7 = 5 ■ 10 - 4 s/d 2 (yellow (light gray) line with tri¬ 
angles). Especially for large external forces, the coarse¬ 
grained rates deviate strongly from their asymptotic val¬ 
ues even for a probe as small as the gold particle. 

The average velocity as shown in Fig. © also strongly 
depends on the friction coefficient of the probe particle, 
especially for large external forces. In this regime, for 
large 7 , the velocity is dominated by the friction experi¬ 
enced by the probe while for small 7 the probe relaxes 
almost immediately and the velocity is dominated by the 
timescale of the motor jumps. 

Another option to reach the fast-bead limit is to use 
very small nucleotide concentrations. In Fig. © we show 
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FIG. 2. (Color online) Coarse-grained rates fl + and fU 
(top) and average velocity (bottom) as functions of / ex d for 
various friction coefficients 7 in the range 5 s/d 2 > 7 > 
5 • 10 - 10 s/d 2 (from bottom to top). With decreasing 7 , 
the rates and the velocity approach the corresponding fast- 
bead limit (solid black lines). Parameters: k = 40d“ 2 , 
c T = c D = 2 ■ 10 - 6 M, cp = 10 - 3 M, Ay = 19, 9+ = 0.1, 
wo exp[/x^?]/c^ q = 3 • 10 7 (Ms)W 


the coarse-grained rates for various ATP and ADP con¬ 
centrations. With decreasing nucleotide concentration 
(at fixed A/i), the rates approach the asymptotic and 
fl~. However, it is very hard to do experiments at con¬ 
centrations smaller than cs 10 -7 M as jumps of the motor 
are then very rare. 

In Fig. [2] and in Fig. [3] the dependence of the coarse¬ 
grained rates on the external force exhibits two different 
regimes. Up to values of the external force of roughly 
15/d, the coarse-grained rates can be well approximated 
by a single exponential dependence on f ex with the same 
slope as in the fast-bead limit, d9 + or dd”, respectively. 
However, for large 7 and large ct, even in this regime, the 
absolute values of the coarse-grained rates deviate up to 
two orders of magnitude from their fast-bead approxima¬ 
tion. For such parameters, assuming a mono-exponential 
dependence on / ex with the above slope would not be 
appropriate either. 

For large external forces, all coarse-grained rates de¬ 
viate significantly from their fast-bead limits. We find 
again a mono-exponential decay for fbut now with 
slope — d whereas grows only linearly with increasing 
/ ex . This so far unaccounted for behaviour can be under¬ 
stood by considering the limit / ex —>• 00 as discussed in 
detail in the Appendix. The crossover from one regime 
to the other occurs beyond the stall force / ex = A y/d. 



f e *d 


f ex d 


FIG. 3. (Color online) Coarse-grained rates fl + and O - as 
functions of / ex d for various ct, cd in the range 2 ■ 10 -5 M > 
Ct, Cd > 2 ■ 1(T 12 M (from top to bottom). With decreasing 
ct,cd, the rates approach the fast-bead limits Cl + and fl _ 
(straight lines). Parameters: k = 40d -2 , 7 = 0.5s/d 2 , cp = 
1(T 3 M, Ay = 19, 9+ = 0.1, w 0 exp[/4 q ]/4 q = 3 ■ 10 7 (Ms) -1 . 


In summary, we find that for the Fi-ATPase under 
realistic experimental conditions the rates in a coarse¬ 
grained description comprising only one effective particle 
that satisfy the LDB condition eq. © and reproduce the 
correct average velocity v can not be written in the form 
of a single exponential dependence on the external force. 


E. Comparison of coarse-grained with averaged 
rates 

Instead of defining the coarse-grained rates according 
to eqs. mm, one might be tempted to use the averaged 
rates @ as a definition for the coarse-grained rates. In 
Fig. we show the averaged rates of our Fi-ATPase 
model as well as their ratio corresponding to the LDB 
condition. We find that both (w + ) and ( w~) (for the lat¬ 
ter less visible in the plot) exhibit non-monotonic depen¬ 
dence on the external force. For external forces slightly 
larger than the stall force, (w + ) increases with increas¬ 
ing / ex due to the fact that in this region the system 
moves backward with motor jumps following the probe 
which leads to a peak at small y in p s ( y ). On the other 
hand, (w~) exhibits a minimum around stall conditions 
for large 7 since in this region, p s (y) misses a peak at 
large y — 1. 

A severe issue appears regarding the LDB condition. 
The corresponding ratio of the averaged rates is also plot¬ 
ted in Fig. [1 where it can be clearly seen that the LDB 
condition is not fulfilled (except in the fast-bead limit). 


F. Without external force 

Even though we have motivated this paper by empha¬ 
sizing that external forces are typically applied to probe 
particles, it should be obvious that our approach holds 
true for molecular motors transporting cargo subject to 
Stokes friction in the absence of external forces. 
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FIG. 4. (Color online) Top: Average rates ( w + ) and (w ) 
as functions of / ex d for various 7 in the range 5 s/d 2 > 7 > 
5 • 1CP 9 s/d 2 . With decreasing 7 , the rates approach f) + , fl~ 
(solid black lines). Bottom: Ratio of + and — rates. In con¬ 
trast to fl + , IP (large red dots), the averaged motor rates do 
not fulfill the LDB condition (solid black line). The parame¬ 
ters are the same as in Fig. [2] 


FIG. 5. (Color online) Coarse-grained rates fl + and IP (top) 
and average velocity (bottom) for various 7 and / ex = 0 
as functions of ct- Since cd and cp are fixed, A/r also 
increases with ct- The rates and the velocity approach 
the fast-bead approximation (solid black lines). Parameters: 
c D = 2 • 10“ 6 M, c p = 1 • 1(T 3 M, k = 40d" 2 , 9+ = 0.1, 
wo exp[/rSj?]/&j? = 3 • 10 7 (Ms) _1 , 7 in the range 5s/d 2 >7 > 
5 • 1CF 9 s/d 2 (from bottom to top). 


For one-particle models, the friction coefficient of the 
probe can not be taken into account explicitly. One 
rather has to incorporate the drag effect of the bead into 
the motor rates [46|. If one wants to analyze experimen¬ 
tal data obtained from probe particles of different sizes, 
one then has to use different values of the motor rates for 
each data set. 

For the rather dilute solutions used in experiments 
[49l |51 . [55'5 one generally assumes that the motor dy¬ 
namics is subject to mass action law kinetics, i.e., that 
the transition rates depend linearly on the correspond¬ 
ing concentration of nucleotides. Obviously, this linear 
dependence holds for all concentrations and beads of all 
sizes for one-particle models. When keeping cd and cp 
fixed, the average velocity of a one-state motor will show 
a purely linear dependence on ct- 

The experimental analysis of the average velocity of the 
Fi-ATPase as function of ct (for fixed cd, cp) reveals a 
saturation of the velocity for large ATP concentrations 
which sets in earlier for large beads [53|. While such a 
saturation is usually attributed to the hydrolysis step, we 
find that a sub-linear dependence of the velocity can also 
be caused by the drag of the probe particle. 

In Fig. 0 the coarse-grained rates as well as the ve¬ 
locity are shown as a function of the ATP concentration. 
With decreasing 7 , the coarse-grained rates approach the 
fast-bead limit and the mass action law kinetics. The ve¬ 


locity is then linear in ct as in a one-particle model. For 
large 7 , eliminating the cargo by coarse-graining yields 
coarse-grained rates that are not linear in the concentra¬ 
tions although the motor rates are still subject to mass 
action law kinetics. Moreover, the velocity then exhibits 
a sub-linear dependence reminiscent of the typical satu¬ 
ration effect for large Ct- 


G. Comparison of full and coarse-grained 
trajectories 

Trajectories of motor and probe generated by a simu¬ 
lation of the complete model of the Fi-ATPase are shown 
in Fig. [G] Additionaly, Fig. [ 6 ] contains a trajectory ob¬ 
tained from simulating the corresponding coarse-grained 
model. The average velocity of both models is the same 
(by definition, see eq. (0), whereas the coarse-grained 
model produces trajectories that are “more random”. 
This behavior occurs since the coarse-grained rates are 
constant (for fixed parameters) and produce a simple bi¬ 
ased random walk. The motor transition rates of the 
complete model, however, depend on the actual position 
of the probe and are therefore implicitly time-dependent. 
Since fast successive motor jumps are suppressed, the tra¬ 
jectory of the complete model is less random |2ll. 1561] . 
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FIG. 6 . (Color online) Trajectories of the one-state model 
for the Fi-ATPase for several parameter sets obtained from 
simulations. The trajectory of the detailed model (motor: 
step-like blue lines, probe: fluctuating red lines) is shown to¬ 
gether with a trajectory of its corresponding coarse-grained 
model (green (light gray)). Parameters: k = 40d -2 , 9 + = 0.1, 
wo exp[/r!p q ]/c^ q = 3 • 10 7 (Ms) -1 , 7 = 0.5s/d 2 , / ex = 0, 
ct = cd = 2 • 10 -6 M, cp = 0.001 M (top left); 7 = 0.5 s/d 2 , 
/ ex = 40 d -1 , c T = c D = 2 • 10 -6 M, c P = 0.001 M (top 
right); 7 = 0.005 s/d 2 , / ex = 0, ct = cd = 2 • 10 -6 M, 
cp = 0.001 M (bottom left); 7 = 0.5 s/d 2 , / ex = 0, ct = 0.001 
M, cd = 2 • 10 -6 M, cp = 0.001 M (bottom right). 


The influence of parameters like the probe size or the 
ATP concentration on the dynamics is visible in the bot¬ 
tom panels of Fig. [6] While the average velocity is almost 
the same, the trajectories of the complete model differ 
significantly. Using a small probe with a small friction 
coefficient, the probe relaxes to the potential minimum 
of the linker before the next motor jump occurs, whereas 
the large probe cannot relax [25} . Large ATP concentra¬ 
tions induce many forward and succesive backward mo¬ 
tor jumps that are absent at lower ATP concentrations. 
These details are not captured in the coarse-grained tra¬ 
jectories. 


III. MOTOR MODELS WITH SEVERAL 
INTERNAL STATES 


A. Explicit motor-bead dynamics and 
coarse-graining procedure 

In this section, we will generalize the model taking into 
account several different internal states of the motor la¬ 
belled by i. The motor states represent the nodes and the 



FIG. 7. (Color online) Network representation of a motor- 
bead model with four internal motor states and discretized 
state space of the probe particle (left). Each row of black dots 
represents one motor state while the dots themselves represent 
specific distances y accessible to the probe particle (via the 
vertical red lines) within the same motor state. Transitions 
between motor states either leave y the same (horizontal green 
lines) or can advance the motor by dfj and change y (diagonal 
blue lines). The top view of this network corresponds to the 
coarse-grained version of this model (right). 


transitions the edges of a network. Transitions between 
the motor states i and j change the free energy by 

A F% = Fj -Ft- A ^ (17) 

where Fj — Fi is the free energy difference of the inter¬ 
nal states of the motor and A/zg- = — A/jT is the free 
energy change of the solvent. Depending on the transi¬ 
tion, A/i" is given by /zt, /zd, Pp or an Y combination 
thereof or 0. Transitions may also advance the motor 
a distance = —d%. Since we allow for several transi- 
tions connecting two states, we assign an additional index 
a to the transitions indicating which link between i and 
j is used. An example for the network of a full system 
comprising motor and probe particle is shown in Fig. [7] 
where the state space of the probe is discretized for better 
presentation. 

The Fokker-Planck-type equation for such models is 
given by 

9tPi(y ) = d y ({d v V(y) - fex)pi(y) + d v Pi{y)) /y 
+ w %(y + d ij)Pj{y + d ij ) - w ij(y)Pi(y)- ( 18 ) 

j,a 

with transition rates of the motor that obey a LDB con¬ 
dition 

jjdSkk-y = expl-acg - V(y + dg) + n»)l. (19) 

The coarse-grained version of such a model should take 
into account the different states of the motor as well as 
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the several possible a-transitions between i and j. Thus, 
the motor network (including all motor cycles) should be 
conserved under coarse-graining. To account for the sev¬ 
eral internal states, we require that the coarse-grained 
rates should obey a LDB condition and the operational 
current [57J from motor state i to motor state j via edge 
a should be conserved. The operational current is the 
sum over all y-dependent net transition currents that 
contribute to the transition i —» j. Conserving the op¬ 
erational currents corresponds to the condition of repro¬ 
ducing the correct mean velocity for the one-state model. 
The above conditions read 


and 


m 


exp[-Ai^ - / ex d“] 


P.O«. _ p.O“ — a 


with the operational current 


( 20 ) 


( 21 ) 


Jll — 


oo 

j [iPi(y)Kj ( y ) - Pj (y + (y + dg)] d y 


= ~3ji 

and the marginal distribution 


( 22 ) 


/ OO 

Pi(y) d y. 

-oo 


(23) 


These equations can be solved for f2g and 0“ : using sim¬ 
ple algebra which yields the rates 


,a = ,a eX P[~A F* - /exdg] 

‘ ij Pie X p[-AF*-f m d%}-P j 

^ = 3?j Pi exp[—— fexdfj] — Pj 


(24) 

(25) 


In principle, it is sufficient to use only eq. (1241) . since 0“ 
takes exactly this form with jg = — jgjAFg = —A 
and dg = —dg. This equivalent procedure would be 
more symmetric and treat all transition rates on an equal 
footing but the LDB condition is then less obvious. Note 
that without the LDB condition (1201) . the stated condi¬ 
tions of Pi and j“- would also be compatible with coarse¬ 
grained rates like the ones in, e.g., [3ll [33j . 

Transitions whose rates are independent of the linker 
elongation y and hence have dg = 0 retrieve their orig¬ 
inal rate constants through this coarse-graining proce¬ 
dure. For such a transition, jg is given by 


3ij = Pi<j ~ p jWji (26) 

with rates fulfilling the LDB condition tug/tug = 
exp[—AFg]. Inserting jg in eqs. (l24l l25l) and using 
the LDB condition and dg = 0 immediately yields 


o“ = w a . 

10 W i-3 ’ 


= w“. 
01- J* 


( 27 ) 


Transitions with rates depending on y but with dg = 0 
have coarse-grained rates that depend on / ex only im- 
plicitely via jg and as will be discussed below in 
section IIII Dl for the chemical transition rates of kinesin. 

The rates determined from the LDB condition eq. (1201) , 
the populations Pi and the operational currents are alge¬ 
braically consistent with the fact that a full set of rates 
fig will uniquely determine the populations Pi on the 
coarse-grained network. Consistency can be seen by in¬ 
tegrating the Fokkcr-Planck equation (fTHl) over y yielding 
the coarse-grained master equation 

d *p> = E iii = E ( 28 ) 

j,a 0 ,ol 

whose stationary solution in the NESS can be expressed 
as a function of the rates fig- [57j, [EU. Thus, the expres¬ 
sion of any current observable in terms of the operational 
currents is consistent with its expression in terms of cycle 
currents on the coarse-grained network. 


B. Time-scale separation 

Similar to the one-state-model, we explore the conse¬ 
quences of a putative time-scale separation between the 
dynamics of motor and probe for each motor transition. 
In the limit 7 — > 0 (formally equivalent to e —> 0 but 
here one would have several within the Fokker-Planck 
equation and all go to 0 ) the solution of eq. (TTH1) in the 
NESS becomes, analogously to [29l. l32|. 

Piiv) = Pi ex p[—P(y) + /exy]/A/\ (29) 

The marginal distribution can be obtained using eq. m 
with its solution for fast bead relaxation 

/ OO 

dtfi(y) dy 

-OO 

= E (E K>» - p «•>») = E E = °- (so) 

j,a j,& 

For Kramers-type transition rates like eqs. CEo] ED, 

Kj ( V) =Ko exp[/x“.’ + -V(y + dg0g’ + ) + V(y)} (31) 
Wji ( y) =kji exp [p?f -V(y- dg0g’ _ ) + V(y)], (32) 

with /zg’ + - /zg’ - = A/zg and fcg/fcg = exp[-Fj + F ,], 
the y-averaged rates (u>g) y and (wfyy become 

«j)y =k?j expiry - /exdg0g’ + ] (33) 

( w ji)y = k ji exp [n°f + /exdgdg’"]• (34) 

The change of chemical free energy A/zg is split into /z“ ,+ 
and /xg’ - indicating that both directions of the transition 
can involve binding and release of the chemical species 
that account for A/zg. The free energy change arising 






from changing the motor state, Fj — Fi, is incorporated 
in the attempt frequencies of the corresponding states. 
Inserting the operational current in the form of eq. m 
with these averaged rates, simple calculus shows that the 
coarse-grained rates (l24l) and (l25l) reduce to 

fig =k?j exp[/x“’ + - /exd“ Cy + ] (35) 

=kfi exp[/x“.+ / ex d“ Off] (36) 

which is again consistent with transition rates of one- 
particle models that assume a purely exponential depen¬ 
dence on the external force. 


C. Example: Fi-ATPase with intermediate step 

1. With external force 

The 120° step of the Fi-ATPase is known to consist of 
a 90° and a 30° substep [53|. Such a stepping behavior 
can be modelled with a unicyclic motor with two internal 
states. A schematic representation of a system compris¬ 
ing a probe particle and a motor with two internal states 
is shown in Fig. [5] The two different pathways for tran¬ 
sitions between the states 1 and 2 correspond to the 90° 
and 30° substeps of the Fi-ATPase, respectively. 

Like in section lll Dl for the one-state model, we examine 
the coarse-grained rates for the 90° and 30° steps and the 
velocity which are shown in Fig. [9] Similar to the 120°- 
scenario, the rates approach their fast-bead limit with 
decreasing 7 . 

As in the one-step model, the dependence of the coarse¬ 
grained rates on the external force shows two regimes. 
For small external forces, the rates can be well approxi¬ 
mated by a single exponential dependence on / ex with 
slope ±d?j6°j in most cases. For large probe parti¬ 
cles, however, the rates neither match the absolute value 
nor show mono-exponential dependence on / ex with the 
above slope. For large forces, the forward rates decay 


«’“(y) u>f?(y + d?i) 



v(t) 


FIG. 8 . (Color online) Schematic representation of a motor- 
bead model for the Fi-ATPase with two inernal states of the 
motor, 1 (blue (dark gray)) and 2 (pale green (light gray)). 
Transition between states 1 to 2 corresponding to the 90° 
(30°) substep are labelled with superscript 90 (30). The tran¬ 
sition rates are chosen accordingly to eqs. (EH EH). 




fexd fexd 



FIG. 9. (Color online) Coarse-grained rates for the 90° (top) 
and the 30° (center) substep and average velocity (bottom) 
as functions of / ex d for various 7 in the range 5s /d 2 > 7 > 
5 • 10~ 10 s/d 2 (from bottom to top). With decreasing 7 , the 
rates and the velocity approach their corresponding fast-bead 
limit (solid black lines). Parameters: k = 40d -2 , ct = cd = 
2 • 1(T 6 M, c p = 10“ 3 M, 0+ i3O = 0.1, exp[/rSj?]/c!j? = 3 • 
10 7 (Ms)-\ fci?ex p^/cd = 3667.5(Ms)- 1 , = 1000s" 1 , 

ki 2 exp[/x^?]/cp q = 40(Msl" 1 . The attempt frequencies are 
chosen on the basis of 15a. 15411 where very small probe particles 
have been used. 


faster whereas the backward rates grow more slowly than 
in the fast-bead limit. 

Concerning the average velocity, strong deviations 
from the fast-bead limit occur only for the largest fric¬ 
tion coefficients. Using small beads, the force-velocity 
relation resulting from our coarse-graining procedure co¬ 
incides well with the one obtained from a one-particle 
model due to the fact that the velocity involves only dif¬ 
ferences of the rates multiplied with the marginal dis¬ 
tribution rather than the rates themselves. For large 
external forces and small 7 , the velocity is significantly 
smaller than in the one-state model since the motor has 
to take two successive steps to cover the full d. The 
force-velocity relations for the two-state as well as for 
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the one-state model reproduce very well the exerimen- 
tally determined force-velocity relation from |[lj for the 
corresponding value of the friction coefficient 7 . 

The limiting cases / ex —> ±00 are more involved here 
than in the one-state model since one has to account 
for the dependence of the P/s on the external force. 
However, as long as the Pj’s do not decay faster than 
exp [—/ ex d“ ], it is still possible to approximate the rates 

mm by 


» -j“ exp[- AF% - fexdfiyPj , (37) 

«“i » -Jij/Pj (38) 


since Pi is bounded by 1 . 

For the Fi-ATPase model, the numerical analysis in 
the / ex —> 00 limit yields a linear dependence of (y) and 
j°j on / ex . We also find that P 2 decays exponentially 
while Pi approaches 1. Hence, and ^21 decay expo¬ 
nentially with slope — di 2 = — 0.75d and — 1 = —0.25d, 
respectively, like in the one-state model but fi®? now 
grows exponentially with a smaller exponent while fl 3 !) 
still grows linearly. 


2. Without external force 

Just as for the one-state-model, we examine the de¬ 
pendence of the coarse-grained rates on the ATP concen¬ 
tration in the absence of external forces. 

Fig. flOl shows the coarse-grained rates for the 90° and 
the 30° substeps as well as the avergage velocity. With 
decreasing 7 , the coarse-grained rates approach the mass 
action law kinetics for the corresponding one-particle 
rates. In contrast to the one-state model, even in this 
limit, the velocity shows saturation. This is due to the 
fact that the timescale of the hydrolysis reaction is in¬ 
dependent of the ATP concentration and represents the 
limiting effect for the velocity. The dependence of the av¬ 
erage velocity on the ATP concentration is reminiscent 
of a Michaelis-Menten kinetics and coincides well with 
experimenal results for several different probe particles 
as shown in (H3 |. 

For large beads, the coarse-graining process yields 
rates that are no longer linear in the corresponding con¬ 
centrations. In this regime, the sub-linear dependence of 
the velocity on the ATP concentration appears already 
for smaller ATP concentrations. Comparing the velocity 
curves of the two-state model with the one-state model, 
we find that for large beads the velocity curves almost 
coincide since in this regime the limiting effect for the 
velocity is the friction experienced by the bead. Thus, 
using large probe particles it is not possible to infer the 
underlying motor dynamics from the characteristics of 
the velocity as a function of the ATP concentration [2f| . 

Fig. E] shows the coarse-grained forward rates for 
three different nucleotide concentrations and for various 
7 chosen as in the experiment [52| . We find that the 90° 




FIG. 10. (Color online) Coarse-grained rates for the 90° (top) 
and the 30° (center) substep and average velocity (bottom) 
for various 7 and / ex = 0 as functions of ct- Since cd and 
cp are fixed, Ay also increases with ct- The rates and the 
velocity approach the fast-bead approximation (solid black 
lines). Parameters: cd = 2 • 10~ 6 M, cp = 1 • 10~ 3 M, k = 
40d -2 , 0 /o, 3 O = 0.1, 7 in the range 5s/d 2 >7 > 5-10 _ 1 °s/d 2 
(from bottom to top). 


rate depends only weakly on 7 for small ATP concen¬ 
trations which is reminiscent of the experimental obser¬ 
vation that the ATP binding rate to the motor depends 
only weakly on the size of the probe (52[. However, for 
large ATP concentrations that were not investigated in 
the experiment, the 90° rate shows a strong dependence 
on 7 . This is due to the fact that for small ATP concen¬ 
trations the relaxation times of all probe particles are in 
the order of, or even faster than, the motor jump rates. 
The results for the 30° rate are consistent with exper¬ 
imental results for the hydrolysis rate [52j| . Increasing 
cp decreases the P; release rate in the experiment as it 
decreases the 30° rate here. 
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FIG. 11. (Color online) Coarse-grained forward rates for var¬ 
ious 7 = 2.1s/d 2 , 0.26s/d 2 , 0.14s/d 2 , 0.017s/d 2 , 0.003s/d 2 , 
0.001s/d 2 , 3.8- 10 _ 4 s/d 2 (from bottom to top) for the param¬ 
eter sets I: ct = 430 nM, cp = 1 nM; II: ct = 1 mM, cp = 1 
nM; III: ct = 1 mM, cp = 200 mM as used in [52|. Parame¬ 
ters that are the same for all sets I-III: cd = 1 nM, k = 40d -2 , 
9g 0 30 = 0.1 and the attempt frequencies fc“- as given in Fig [9] 
The values of cd = cp = 1 nM are a rough estimate because 
there is no information about these concentrations in [52] . 


D. Example: Kinesin 



FIG. 12. (Color online) 6 -state-model representing a kinesin 
motor adapted from [jj. The transition between states 2 and 
5 is purely mechanical and corresponds to a step of length 
d whereas all other transitions are pure chemical transitions. 
The motor model includes three cycles: T which, in the + 
direction, includes ATP hydrolysis and forward stepping, B 
which includes ATP hydrolysis and backward stepping in its 
+ direction and a pure chemical cycle (around the circle) that 
includes hydrolysis/synthesis of two ATP. 


As a final more complex example, we apply our coarse- 
graining method to a model with a multi-state motor. We 
choose the well-studied 6 -state-model representing a ki¬ 
nesin motor introduced in jBJ, see Fig. 1121 Implementing 
the probe particle and an elastic linker V(y), we adopt 
the transition rates of the motor from [Bj and replace the 
dependence on the external force by the dependence on 
the elongation of the linker, 


w % 5 (v) = fc 25 exp [~V(y + dd + ) + V{y)} : 

W 52 (v) = k 52 exp[—y (y - dd~) +V(y)\, 

+ 2 exp [ntj] 

«^,chem - ii 1 + exp [d y V{y) X ij] ’ 


2exp[M y ] 

,chem *l + exp [dyV (y) Xij ]' 


(39) 

(40) 

(41) 

(42) 


The hrst two rates belong to the mechanical transition, 
the lower two rates represent the chemical transitions 
which depend on the instantaneous force exerted by the 
linker with a chemical load-sharing factor Xij • see [S|- 
The change of chemical free energy yf^ = y ^, y ^, yp de¬ 
pends on which transition involves binding of the corre¬ 
sponding nucleotide. We choose again V(y) = ny 2 / 2. 

The coarse-grained rates for the mechanical transition 
are shown in Fig. [TIJ] With decreasing 7 , the rates 
approach their fast-bead limit which corresponds to the 
rates used in [5| while strong deviations occur for finite 7 
especially for assisting external forces. The friction coeffi¬ 
cient of a probe of size 500 nm as in Q can be calculated 
using Stokes’ law yielding 7 ~ 7.7 • 10 - 5 d 2 /s. For fric¬ 
tion coefficients in this range (light green (light gray) line 
with triangles), our coarse-grained rates show a distinct 
deviation from the one-particle rates (solid black lines). 



ex 14 ,/ex 1 


FIG. 13. (Color online) Coarse-grained rates for the me¬ 
chanical transitions (with y-dependence) for various 7 in 
the range 0.077s/d 2 > 7 > 7.7 • 10~ lo s/d 2 (from bot¬ 
tom to top). The rates approach the one-particle rates 
from [B] (solid black lines). Parameters: k = 10d~ 2 [BBI ]. 
ct = 0.001M, cd = cp = 10 _9 M (estimated), 9 + = 
0.65, Xij = 0.25,0.15, fci2 exp[y^?]/cT q = k4,5exp[y^]/c^ = 

2 ■ 10® (Ms) 1 , &21 = &23 = &34 = &56 = &61 = 
100(s) _i , fc 32 exp[^]/cD q = fees exp[/X^]/Cn q = 2 • 10 4 (Ms) _1 , 
fc 4 3 exp[/ip q ]/cp q = fei6 exp[yp q ]/cp q = 2 • 10 4 (Ms)^ 1 , /C 25 = 

3 • 10 5 (s) 1 , k $2 = 0.24(s) 4 , &54 = (^ 52 /fos ) 2 ^21 - 


However, the average velocity (obtained from our coarse¬ 
grained rates) as function of the external force coincides 
very well for almost all 7 with the velocity curve obtained 
from the bare motor model, see Fig. [14] Like for the 
Fi-ATPase model discussed in section UlI Cl this agree¬ 
ment is due to the fact that the velocity involves only the 
difference of the rates multiplied with the marginal dis¬ 
tribution. If one investigates only force-velocity curves, 
the discrepancies between the coarse-grained rates and 
the one-particle rates are hardly visible. 

In contrast to the coarse-grained rates of the Fi- 
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FIG. 14. (Color online) Left: Average velocity for the model 
with probe particle (colored lines with symbols) an the one- 
particle model from J5| (solid black line). Right: Coarse¬ 
grained rates for chemical transitions (with y-dependence) for 
7 = 0.077s/d 2 . Other parameters as given in Fig. 1131 



ATPase models, the coarse-grained rates for the mechan¬ 
ical transition of the kinesin model show more struc¬ 
ture especially for negative, i.e., assisting external forces. 
Since the kinesin model contains several internal motor 
cycles, depending on the external force the dominant cy¬ 
cle can change leading to crossover regimes with changing 
weight of the probabilities Pi. 

The dependence of the coarse-grained rates for chem¬ 
ical transitions on the external force is visible in Fig. 
M Although there is no explicit dependence on external 
forces for pure chemical rates since =0, /“■ and Pi de¬ 
pend on / ex via y. The operational current for transitions 
within the .F-cycle in + direction decreases with increas¬ 
ing / ex whereas the operational currents within the £>- 
cycle in + direction slightly increase with / ex which can 
be explained intuitively since the motor prefers “back¬ 
ward” cycles for large opposing forces. However, all 
coarse-grained rates decrease with increasing / ex simi¬ 
lar to the bare motor rates (Idll l42l) which decrease with 
larger y ; a situation that is more likely to appear for large 
external forces. 


IV. EXPERIMENTAL IMPLEMENTATION 

In order to practically apply the coarse-grained de¬ 
scription, one has to determine the marginal distribu¬ 
tions Pi, the operational currents _)“■ and the free en¬ 
ergy differences A F^. For multi-state motors, this is a 
rather challenging task since only a few quantities can 
be extracted reliably from the experimentally measured 
trajectory of the probe. Note, however, that this prob¬ 
lem does not happen exclusively in our approach but is 
inevitable whatever method is used to infer motor prop¬ 
erties from such trajectories. 

In the following, using the 90-30 model for the Fi- 
ATPase, we will illustrate how these quantities can be 
estimated. If all motor transitions involve mechanical 
transitions with different step sizes, the plateaus in the 
probe trajectory can be assigned to specific correspond¬ 


ing motor states. Since after a large enough time inter¬ 
val all possible transitions will have occured, one is also 
able to reconstruct the links connecting the states. The 
marginal distributions Pi are then given as the fraction of 
time that the corresponding motor state is occupied. For 
the two-state model of the Fi-ATPase, we assign plateaus 
in the probe trajectory that are followed by a fast 90° for¬ 
ward or 30° backward displacement to motor state i = 1 
and plateaus that are followed by fast 90° backward or 
30° forward displacement to i = 2. In principle, there 
are several possibilities to reconstruct hidden variables 
from partially visible trajectories [59l46l| . Here, we will 
use a simple algorithm which sets i = 2 if four consec¬ 
utive data points are within a specific range around 90° 
and otherwise i = 1. The marginal distributions Pi, P 2 
are then represented by the fraction of data points with 
assigned i = 1,2. 

If the motor is not very complex, the operational cur¬ 
rents /“■ can be obtained rather easily since they are pre¬ 
cisely the net currents between two motor states. For 
unicyclic motors, all operational currents are equal to the 
average velocity divided by d, the operational current of 
an ATP binding transition is the net disappearance rate 
of ATP in the solution (given that there are no other ATP 
binding reactions), etc.. If all motor transitions involve 
mechanical transitions with different step sizes, the oper¬ 
ational currents between any two states can be obtained 
by counting the number of transitions of a specific step 
size from i —> j, nfj, and j —> i, n/j. The (time) average 
of this current using one long trajectory of length f to t is 
then given by 


iij=( n ij- n ji)/ t t ot- (43) 

In our example, in order to estimate jf$ we have to count 
the number of sudden displacements of “size” 90° either 
from the trajectory of the probe directly or from the re¬ 
constructed trajectory of the motor using the assignment 
rule mentioned above. If the time resolution of the tra¬ 
jectory is very coarse or if the reconstruction method is 
rather inaccurate, jumps that consist of fast consecutive 
90 and 30 jumps with apparent step size 120° will appear 
which have to be included in the number of 90° (and also 
30°) jumps. Fig. fl5l shows a reconstructed motor trajec¬ 
tory obtained with the algorithm mentioned above. We 
have used a trajectory of the probe from our simulations 
as “experimental data”. Compared to the original motor 
trajectory, this reconstruction captures the average dy¬ 
namics quite well. Large fluctuations of the probe can 
generate additional apparent motor jumps in the recon¬ 
structed trajectory that are absent in the original one. 

Finally, the estimation of the free energy difference 
A F°- = Fj - Fi — A/r“- is slightly more involved. In 
equilibrium (Ay, = 0, / ex = 0), detailed balance holds, 

= p?(y + d ?j) (u] 

™%(y + d W p?{y) ’ 1 j 

with the Boltzmann distribution p® q (y) = 
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P® q exp[— V(y)\/Af. Inserting this expression yields 


Pp/Pr = exp [-Fj + Ft + A/r“’ eq ] = exp[-AP“] 

(45) 


for the marginal distributions in equilibrium. Note that 
A/r“-’ eq 0 if the corresponding transition comprises only 
binding or release of nucleotides. Thus, the equilibrium 
free energy difference A(which explicitely depends 
on the equilibrium concentrations) can be obtained from 
the ratio of the marginal distributions under equilibrium 
conditions. Using /x* = /x® q + ln(cj/c® q ), we find that 

AF“ = AJ-“±^ln^ (46) 

k 


with k = T, D, P and the sign depending on which bind¬ 
ing or release event corresponds to the transition ij, a 
[62]. Hence, the free energy difference A F-j needed for 
the coarse-grained rates can be expressed by the equi¬ 
librium free energy difference AJ-T obtained from ex¬ 
perimental data at equilibrium conditions and the nu¬ 
cleotide concentrations with respect to the equilibrium 
concentrations corresponding to the conditions used to 
obtain AJq“. For the 90-30 model, we have P| q /P® q = 
expt-APfO] = exp[—APf))] with -APfo = -P 2 + Pi + 
— fp = —P 2 + Pi + = —A Fp. since Afi = 0 in 

equilibrium. 

Once these quantities have been estimated, there are 
no additional fit parameters needed or left. All concen¬ 
trations as well as the external force are usually known 
from the experimental setup. To obtain the coarse¬ 
grained rates from the probe trajectory of our 90-30 
model, we then proceed as follows. First, we choose 
equilibrium conditions and obtain APi 2 from the ra¬ 
tio of marginal distributions. Then, we change to non¬ 
equilibrium concentrations and estimate Pi, P 2 and the 
operational current jf®- The coarse-grained rates, ac¬ 
cording to eqs. (l24l 1251) . are then given by 


o90 -90 c T CD q exp[-APi 2 -/ ex d^]/(4 q c D ) 

12 Jl2 Pic T c e D q exp[-APi 2 - / ex d?°]/(4 q c D ) - P 2 ’ 

(47) 

q90 _ -9C_ l _ 

21 Jl2 PiCtcS 1 exp[-APi 2 - / ex d?°]/(4 q c D ) - P 2 ’ 

(48) 

o 30 _ „-90 c P q exp[APi 2 — / ex d 2 i]/cp . . 

21 12 P 2 Cp q exp[AP 12 - / ex d|?]/c P - Pi ’ 

^12 = Jl2 75 eq 7 ^30l / 7T' (^0) 

P 2 Cp q exp]APi 2 - /e X d 2 ij/cp - Pi 

A comparison of the coarse-grained rates and related 
quantities obtained from the full theoretical model and 
from the reconstructed one estimated using the probe 
trajectory is shown in table Q] We find quite good agree¬ 
ment between the original and the reconstructed quan¬ 
tities with a maximum error of 14% except for the Q'P 
rates which have a maximum error of 24%. 



FIG. 15. (Color online) Comparison of the simulated trajec¬ 
tory of motor and probe with a trajectory of the probe and 
the estimated motor position that was reconstructed using the 
simulated trajectory of the probe. The trajectories are shifted 
for better visibility. Parameters: k = 40 dT 2 , 7 = 0.005s/d 2 , 
c T = c D = 2 • 10 _ 6 M, c p = 10 - 3 M, &r q = 3.33 ■ 10 _ 7 M, 
Cq 1 = 0.0682M, Cp q = 1M, / ex = 0, lower boundary to set 
i — 2: x — L*J = 0.375d, upper boundary to set i = 2: 
x — [x\ = 0.89d. 


TABLE I. Comparison of the coarse-grained rates and other 
relevant quantities obtained from the simulation of the full 
model with the ones estimated using the reconstructed mo¬ 
tor trajectory. The trajectory used to obtain these values is 
shown in Fig. [15] 



simulation 

estimate 

Pi 

0.944 

0.952 

Pi 

0.056 

0.048 

jll [1/S] 

52.292 

52.246 

AJka 

3.216 

3.140 

ni° 2 [1/s] 

55.325 

54.899 

nlS [i/s] 

0.00676 

0.00620 

U 21 [l/s] 

937.1 

1082.37 

U 12 [1/s] 

0.037 

0.046 


The 90-30 model thus provides a useful demonstration 
of the experimental applicability of the coarse-graining 
method showing that it is possible to estimate the coarse¬ 
grained rates from experimental accessible data if the un¬ 
derlying motor network is not too complex. Considering 
the simplicity of the applied reconstruction method, the 
accuracy of the estimates is rather encouraging. 


V. INVARIANCE OF ENTROPY PRODUCTION 
AND EFFICIENCY 

An important question for any coarse-graining method 
concerns its effect on entropy production. In general, 
a coarse-grained description without imposed time-scale 
separation or detailed balance for the eliminated vari¬ 
ables often underestimates the entropy production of the 
system [2914321 . l40( | . In this section, we show that for 
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the type of models considered here, our coarse-graining 
method conserves the entropy production even if there 
is no time-scale separation between the eliminated and 
remaining degree of freedom. 

Since transitions can be uniquely attributed to motor 
or probe particle, the total entropy production of the sys¬ 
tem [28| can be split in two parts, analogous to bipartite 
or partially masked systems |63|, |64| , 


= 


oo 

■£/ 


r 7 jf{y) 

J- oo Pi(y) 

Pi(y)Kj(y ) ln 


d y 


Pi{y)wii(y) 


pj{y + d ii) w u(y + d ?j) 


= si 


tot 


d y 


(51) 


entropy production (1531) plus a contribution of the mi¬ 
crostates corresponding to a mesostate (which are elimi¬ 
nated during coarse-graining) plus a contribution due to 
the fact that jumps between mesostates can occur involv¬ 
ing different microstates. In our framework, the total 
entropy production is already recovered by the coarse¬ 
grained entropy production. The two additional contri¬ 
butions which correspond to the total entropy production 
of the probe particle and the average total entropy pro¬ 
duction of the motor minus the coarse-grained entropy 
production cancel each other. 

We finally show that our coarse-graining procedure 
also preserves the energy transduction, or thermody¬ 
namic, efficiency rjr defined as the ratio of the extractable 
jjower Wout and the rate of chemical energy input A y 


where jf(y ) = {{d v V(y) - fex)pi(y ) + d y pi(y ))/7 is the 
current due to the motion of only the bead for fixed i. 
Obviously, both S^ ot and 5 1 ™,. are non-negative. 

The total entropy production (1511) can be calculated 
using the LDB condition (fTHl) as 


/ OO 

(d y V(y) - fex)ji(y)dy 

t -°° 

/ OO 

Pi(y)Kj(y ) i A P-ij - F o + p i 

-00 

-V(y + d°j) + V(y))dy 
= E - /ex» > 0. (52) 


i<j,a 


Using partial integration, it can be easily seen that the 
parts involving V(y) cancel, i.e., the energy of the linker 
is constant on average. The total entropy production is 
then given by the chemical free energy consumption that 
is not transformed into mechanical power. 

For the coarse-grained description, the total entropy 
production contains only contributions from the effective 
jump process, 


^to S t = E 111 S’ = E ln (53) 

i,j,a J* 




Using the LDB condition for the coarse-grained rates (®D 
and the condition on the operational current (1211) yields 


s t c 0 s t = e - /ew (54) 

i<j,QC 

which is precisely (1571) . For these models for which the 
state space of the eliminated degree of freedom does not 
contain entropy producing internal cycles, the average 
total entropy production in the NESS remains invariant 
under our coarse-graining procedure. 

It is also instructive to apply the entropy-splitting 
scheme introduced in |3lj ] to our coarse-graining proce¬ 
dure. In [31], it was shown that the total entropy pro¬ 
duction can be written as a sum of the coarse-grained 


Vt = 


Wout 

Ay 


(55) 


For the systems we have studied so far, as long as the 
external force is smaller than the stall force, the power 
output is given by W out = f ex v and the power input by 
J2i<j a A dijJij which leads to the efficiency 


Vt 


fexV 




A PijJ 


Oi 

ij 


(56) 


which is the same in the coarse-grained description since 
v, j°j and A/z“ are conserved. 

For motor models with tight coupling or multi-state 
models with a single motor cycle, the rate of chemical 
energy input equals the velocity Ay, = vAy/d and the 
efficiency reduces to t]t = /ex/ Ay. In general, however, 
any idle cycles of the motor increases the rate of chemi¬ 
cal input over the velocity and therefore reduces the effi¬ 
ciency. 


VI. STALL FORCE AND RATE ANOMALY 

Coarse-graining multicyclic motor models as developed 
here reveals a remarkable feature concerning the stall 
force with significant implications on the interpretation 
of experimental data. For an example, consider the ki- 
nesin motor for the parameters chosen in Fig. 1131 Fig. 
m shows that the stall force is a function of the size of 
the attached probe particle. Generally speaking, the stall 
force can indeed depend on the size of the probe since the 
network of the full system comprises more cycles than the 
coarse-grained or bare motor network, see Fig. [7] Vary¬ 
ing the size of the probe, the relative weight of the cycles 
in the full system and hence their contribution to an oper¬ 
ational current can change yielding a varying stall force. 
Thus, the experimentally obtained stall force corresponds 
to the stall conditions of the motor-probe complex but 
does not necessarily represent the stall conditions of the 
bare motor. If one is interested in the latter one should 
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FIG. 16. (Color online) Average velocity of the kinesin model 
as a function of the external force. Depending on the size of 
the probe, stall conditions are reached for different values of 
the external force. Parameters as in Fig. 1131 

use very small probe particles since the limit of vanishing 
friction coefficient 7 is equivalent to applying the force di¬ 
rectly to the motor. As discussed below, the stall force 
is independent of 7 for one-state or unicyclic multi-state 
motor models. Hence, an experimentally observed varia¬ 
tion of the stall force with probe size can be used as proof 
that the motor is indeed multicyclic. 

The varying stall force has also implications on the 
transition rates. In Fig. [13] a close look around f ex d = 14 
shows that these data points are missing for the follow¬ 
ing reason. For all investigated models, we find that 
if, as a function of the external force, the sign change 
of an operational current depends on the friction coeffi¬ 
cient 7 , the coarse-grained rates corresponding to this 
transitions can become piecewise negative. This phe¬ 
nomenon occurs when the affinity of the affected tran¬ 
sitions, In /(Pjtl'ji )], has the opposite sign of jJT 
An isolated sign change in the denominator of eqs. (IM1 
E5l) leads to a pole in the corresponding rate. Such an 
anomaly in f 1 ^- necessarily implies a corresponding one 
in since the ratio of the effective rates obeys the lo¬ 
cal detailed balance condition which enforces the same 
sign for both rates. In this range, the coarse-graining 
scheme introduced here fails to produce physically ac¬ 
ceptable rates. In practice, one should discard the results 
at least when either a rate is negative or becomes larger 
than the rate for vanishing bead size. In Fig. [TT] where 
we zoom into the range around the stall force, this range 
is shaded gray. Taken at face value, this phenomenon 
looks like a short-coming of our approach. It is the price 
to pay for requiring over the full parameter range both 
the local detailed balance condition and the correct net 
currents from any one motor state to any other. While 
the negative rates do not allow for a sensible physical 
interpretation, they can nevertheless be used to calcu¬ 
late average quantities and yield, e.g., the correct entropy 
production as shown in section El 

This stall force anomaly with a corresponding range 
of negative rates occurs neither for any one-state motor 
model nor for unicyclic motors around stall conditions 
since only one motor cycle contributes to all cycles of the 



/exd /ex d 

FIG. 17. (Color online) Detail of the coarse-grained mechan¬ 
ical transition rates of the kinesin model as shown in Fig. 1131 
Near the stall force at / ex d — 14 these rates exhibit a pole. 
In the gray shaded range, they should not be interpreted as 
physical transition rates. 

full system causing the zero of and In [Pi^lfj / (PjCl^)] 
to occur for the same / ex . We also found several mul¬ 
ticyclic motors that do not lead to negative rates, e.g., 
the kinesin model if one assumes the chemical rates to 
be independent of y. A derivation of the precise condi¬ 
tions under which for multi-cyclic motor models a pair 
of effective rates diverges or becomes negative must be 
left to future work. We stress, however, that in all ex¬ 
amples shown in this study, this anomaly occurs only in 
the narrow range shown in Fig. 1171 From a practical 
point of view, it may therefore not be as relevant as it is 
intriguing from a theoretical perspective. 


VII. CONCLUSION 

Most experiments on molecular motors comprise some 
kind of probe particle. Therefore, any theoretical mod¬ 
elling with parameters estimated from experimental data 
will explicitly or implicitly contain characteristics of the 
probe particle. 

In this paper, we have introduced a systematic coarse- 
graining method that allows to reduce motor-bead mod¬ 
els to effective one-particle motor models. This coarse- 
graining procedure provides a compromise between a one- 
particle description that is simple to handle and a de¬ 
tailed model comprising the dynamics of the full sys¬ 
tem. It yields an effective one-particle model maintain¬ 
ing the true motor network, where the influence of the 
probe is naturally incorporated without any additional 
assumptions since the simplification of the description 
takes place a posteriori. Any external force acting on the 
probe is then acting on the effective motor directly. The 
coarse-grained rates obey a LDB condition and yield the 
correct net currents. Fixing the marginal distribution 
and the average currents, there is still freedom on how 
to choose the rates. Only with the LDB condition the 
effective rates are determined uniquely. 

Applying the coarse-graining procedure to motor-bead 
models, we find that in general the coarse-grained rates 
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do not show a single exponential dependence on the ex¬ 
ternal force in contrast to what is often assumed for me¬ 
chanical transition rates in one-particle models. Only in 
the often unrealistic limit of fast bead relaxation, the 
coarse-grained rates reduce to the corresponding one- 
particle rates. 

In the absence of external forces, in general the coarse¬ 
grained rates are not proportional to the ATP concen¬ 
tration even if the motor rates obey mass action law ki¬ 
netics. This feature originates from the drag effect of 
probe (due to friction) that is incorporated in the coarse¬ 
grained rates. For the same reason, the average velocity 
shows a sub-linear dependence on the ATP concentration 
even for a one-state motor model. Assuming an a priori 
one-state model with external force acting directly on the 
motor, one would have either to use a rather counterin¬ 
tuitive complex force-dependence of the transition rates 
or to introduce additional motor states in order to obtain 
a sublinearly growing velocity caused by the drag of the 
probe. In a one-particle description, the effect of large 
probe particles on the dwell time distributions could also 
be mistaken as signature of additional motor states thus 
leading to an overly complex motor network (2lL [56} . 

Considering the influence of the coarse-graining proce¬ 
dure on the stochastic thermodynamics of the system, we 
show that the total entropy production remains invariant 
under coarse-graining. This is due to the fact that, on 
the one hand, the state space of the eliminated degree of 
freedom contains no entropy producing cycles. On the 
other hand the design of the coarse-graining procedure 
is also important. It has to conserve the motor network 
as well as the net currents and provide transition rates 
fulfilling a LDB condition. Likewise, the thermodynamic 
efficiency remains invariant in our scheme. 

Our coarse-graining method conserves average quanti¬ 
ties like the entropy production or operational currents 
although eliminating the dynamics of the probe particle 
strongly affects the cycle structure of the full system. In 
order to preserve also fluctuations of current observables 
in the long-time limit it was found that coarse-graining 
methods should conserve the cycle structure of the full 
system pffil. 137} . 

From the experimental point of view, in order to obtain 
the simpler effective model, the underlying mesoscopic 
modelling need not to be known since all these quanti¬ 
ties enter the coarse-grained description via the net cur¬ 
rents and the marginal distributions which, in principle, 
can be extracted from the experimental data as we have 
demonstrated using a two step model for the Fi-ATPase. 

The main advantage of the coarse-graining procedure 
introduced here is that once the rates have been obtained 
from experimentally accessible quantities, they automati¬ 
cally fulfill a LDB condition and provide the correct aver¬ 
age currents, i.e., velocity, entropy production, hydrolysis 
rate etc.. 

For multicyclic motors, the coarse-graining procedure 
can yield rates that can have poles and become (piece- 
wise) negative. If this scenario occurs, the coarse-grained 


rates lack a physical interpretation as transition proba¬ 
bilities in this range but they can still be used to calculate 
average quantities. For this class of motors the stall force 
typically depends on the size of the probe particle, i.e., 
the friction coefficient. Applying naively a one-particle 
model to such an experimental setup would not allow 
to determine the energy transduction mechanism of the 
motor correctly. For one-state motors, the coarse-grained 
rates are always positive. 

So far, we have discussed coarse-graining only under 
NESS conditions. In principle, the coarse-graining pro¬ 
cedure as introduced in sections III Bl and IIII Al can also 
be applied to non-stationary states, e.g., if the nucleotide 
concentrations are not constant and A/t decreases with 
time @,@|. Such a scenario would yield time-dependent 
Pi’s, net currents, LDB conditions and therefore also 
time-dependent coarse-grained rates. 

Further generalizations might include other types of 
models representing the full system. While developed 
here for discrete motor models, the coarse-graining pro¬ 
cedure should be also applicable to continuous motors 
moving in a tilted periodic potential where the poten¬ 
tial minima will become the discrete states of the coarse¬ 
grained effective motor. The introduction of the index a 
in principle also accounts for more involved potentials or 
free energy surfaces that depend on both the motor and 
the probe state. 
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Appendix: Limiting case: Large applied force 

In the limit of large external forces, / ex —>• 00 , the 
coarse-grained rates mm can be expressed as 

ss — vexp[Ay — f ex d\/d (A.l) 

k. — v/d. (A.2) 

While A/i is independent of the external force, the aver¬ 
age velocity is a function of the external force 

v = ( d v V ( y) - / ex )/7 = k( 2/}/7 - /ex/ 7 - (A.3) 

It becomes negative for forces larger than the stall force 
A/i/d which ensures that both f 1 + and f 1~ are positive. 
If there is no time-scale separation between the dynamics 
of motor and probe, (y) grows linearly in / ex for / ex —>• 00 
with a smaller slope than 1/k. On the other hand with 
time scale separation, we have ( y) = / ex /fc. Note that 
within time scale separation, the average velocity has 
to be calculated using the average velocity of the mo¬ 
tor, eq. 0 : since the “average velocity” of the probe 
(d y V(y) — / ex )/7 is zero as a result of the fast-bead limit 
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of eq. (0. Due to the linear dependence of (y) on f ex , 
the average velocity, and therefore also f l~, are then pro¬ 
portional to the external force whereas the exponential 
factor dominates for Q + , 


~ /ex exp[A/r - f ex d]/ (yd) (A.4) 

~ fex/ild). (A.5) 

In the opposite limit of a large assisting force / ex —> 
—oo, the coarse-grained rates 0® become 

Q+ ttv/d (A.6) 

fl - «vexp[—A n + f ex d\/d (A.7) 


As above, the average y grows linearly and the velocity 
is proportional to f ex if there is no time-scale separation 
which leads to 

~ |/ex|/(7d) (A.8) 

~ I/ex | exp[—A/x + /exd]/(yd). (A.9) 

This simple analysis clearly shows that the coarse-grained 
rates do not coincide with the often a priori assumed 
single exponential force-dependence of one-particle rates. 
Within our numerical analysis, the asymptotic behaviour 
appears for |/ ex | > 500/d. The regime for large forces 
shown in Figs. [2] and [3] is not the asymptotics yet. How¬ 
ever, since ( y ) is also linear in f ex in this region yet with 
different slope, v is still proportional to f ex . 
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